NACA TN 3467 



NATIONAL ADVISORY COMMITTEE "g 
FOR AERONAUTICS ^ 

TECHNICAL NOTE 3467 


EFFECT OF INTERACTION ON LANDING -GEAR BEHAVIOR AND 

DYNAMIC LOADS IN A FLEXIBLE AIRPLANE STRUCTURE 

By Francis E. Cook and Benjamin Milwitzky 

Langley Aeronautical Laboratory 
Langley Field, Va. 


NACA 


Washington 
August 1955 


T 

(. 

j 


mm 

TECKJ3';GA1. ! "■ 


r>* ^ 

t‘ u, 


TECH LIBRARY KAPB, NM 


NATIONAL ADVISORY CCMilTUEE FOR AERONAUTICS 


TECH UBRARY KAFB, NM 



QDbbSVM 


TECHNICAL NOTE 5^67 


EFFECT OF INTERACDION ON LANDING-SEAR BEHAVIOR AND 
DYNAMIC LOADS IN A FLEXIBLE AIRELANE STRUCTURE 
By Francis E. Cook and Benjamin Milwitzky 


STMMARY 


The effects of interaction between a landing gear and a flexible 
airplane structure on the behavior of the landing gear and the loads in 
the structure have been studied by treating the equations of motion of 
the airplane and the landing gear as a coupled system. The landing gear 
is considered to have nonlinear characteristics typical of conventional 
gears j namely, velocity-squared damping, polytropic air-compression 
springing, and exponential tire force-deflection characteristics. For 
the case where only two modes of the structure are considered, an equiv- 
alent three-mass system is derived for representing the airpl^e and 
landing-gear combination, which may be used to simulate the effects of 
structural flexibility in jig drop tests of landing gears. 

As examples to illustrate the effects of interaction, numerical cal- 
culations, based on the structural properties of two large airplanes 
having considerably different mass and flexibility characteristics, are 
presented. For the particular caaes considered, it was found that the 
effects of interaction can result in appreciable reductions in the magni- 
tiide of the landing-gear force, particularly -when the flexibility of the 
airplane structure is large and the natiaral frequency is small. Thus, 
neglect of interaction effects, that is, the use of the landing-gear 
forcing function for a rigid airplane, in a dynamic analysis of a flexible 
airplane can lead to the calcxilation of excessive loads in the airplane 
structure . 

In the case of one of the airplanes considered, che structiural loads 
calculated from the interaction solutions are ^zeater than those for a 
cotopletely rigid airplane treatment (rig'’'^ structure subjected to rigid- 
body forcing function) because the effects of dynamic magnification more 
than overcome the reduction in landing-gear force due to interaction. 

In the case of the second airplane, br a-ise of the relatively large natu- 
ral period of the structure in comparison with the duration of the impact 
pulse, the dynamic magnification factor is appreciably less than unity. 
This effect, coupled with the reductions in landing-gear force due to 
interaction, results in structural loads that are less than those for a 
rigid airplane. 
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INTRODUCTION 


In the design of landing gears It is usually assumed that the air- 
plane is a rigid body and development tests are frequently carried out 
in a drop-test Jig with a landing gear attached to a concentrated mass. 

In so doing, it is tacitly assumed that the interaction between the 
motions of the landing gear and the deformtlons of the airplane structure 
has little or no effect on the behavior of the landing gear. Also, load 
time histories obtained on a rigid-body basis are often used as the 
forcing function in a dynamic analysis to determine the inertia loads 
and stresses in flexible airplane structures, again under the assumption 
that the behavior of the landing gear is Independent of the effects of 
airplane flexibility. Although it has been recognized that this assianp- 
tion is not altogether valid, the errors involved have not been considered 
particularly significant in the past because: (a) The errors were thought 

to be on the conservative side and (b) until comparatively recently main 
landing gears have generally been located very close to the nodes of the 
fundamental bending mode of the wing, and the airplane therefore closely 
approximated a rigid body insofar as the behavior of the landing gear is 
concerned. However, the trend toward increased size of airplanes, the 
disposition of large concentrated masses in outboard locations in the 
wings, the use of thinner wings, and the development of unconventional 
configurations tend to increase the flexibility of the airplane structure 
and reduce the natural frequencies of vibration. These characteristics 
tend to cause an Increase in the amplitudes of the oscillatory motions of 
the landing-gear attachment points relative to the center of gravity of 
the flexible system diurlng intact so that the effects of interaction are 
Increased, both with regard to the behavior of the landing gear and the 
dynamic loads in the structure, particularly when the natural period of 
the fxindamental mode of the structure approaches the time duration of 
the lopact pulse. 

A number of analytical sttidles and some simplified model tests 
(refs. 1 to 5) which have been made to evaluate the effects of structiural 
flexibility on landing-gear loads have indicated some reduction in landing - 
gear force due to the effects of structiiral deformation. However, in 
view of the fact that these previoias Investigations considered only rather 
highly Idealized linear-spring landing gears with either no damping at 
all or viscous dandling, a further study of the effects of interaction 
between the landing gear and the airplane structure has been made with 
a more realistic representation of the landing gear. In the present 
analysis, as in reference 6 , the landing gear is considered to have 
velocity-squared damping, polytropic air-compression springing, and 
exponential tire force-deflection characteristics, as is the case with 
conventional oleo-pneimiatlo landing gears in current use. The particular 
purposes of this investigation are to evaluate the effects of interaction 
on landing-gear behavior and to study the errors introduced into the cal- 
culated loads in the structure (applied loads, accelerations, bending 
moments, and shears) when a dynamic analysis is made on the basis of 
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applying the landing-gear forcing function for a rigid body to a flexible 
airplane. For these purposes, case-history studies, based on the struc- 
tiiral properties of two large airplanes having considerably different mass 
and flexibility characteristics, are presented. In order to cover a range 
of parameters, the landing gear of each airplane was assumed to be located 
at three arbitrary spanwise positions in addition to its original location. 

The basic analysis of the landing gear and the airplane structure 
as a coupled system is presented in a general form. In the numerical 
examples presented, however, the system is simplified by considering the 
motions of the airplane in its first two structural modes only. With 
these restrictions, the combination of airplane and landing gear can also 
be represented by an equivalent three-mass system which may be used in 
Jig drop tests of landing gears to siniulate the primary effects of struc- 
tural flexibility. A similar type of concentrated-mass system was used 
in the study of the hydrodynamic impact of a flexible seaplane in ref- 
erence 7* 


SYMBOLS 

General 


g gravitational constant 

t time after initial contact 

T time variable of Integration 

T time to ma-iriTmm landing-gear force 

t’ time after maximum landing-gear force 

VVq vertical velocity at Initial contact 

ft circular frequency of sine pulse 

ftq circular frequency of cosine pulse 

A any variable 

value of any variable A at end of pth Interval subsequent 
to beginning of nianerical-integration procedure 
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Aa 

Ah 

Al 

An 

Aq 

Ap 

Cd 

^Vb 

^g(^) 

“u 

m,r,ni’ 

n 

Pao 

P 

s 

•fcl 

0 

Wu 

Zf 


Landing Gear 
pnenniatic area of shock stmt 
hydraulic area of shock strut, ,Ai^ - Ag 

* I-* 

internal cross-sectional area of ' shock-strut Inner cylinder 

net orifice area of shock strut, Aq - Ag 

area of fixed opening in orifice plate 

cross-sections, 1 area of metering pin or rod in plane of 
orifice 

orifice discharge coefficient 

vertical component of force in shock strut subsequent to 
beginning of shock-strvct deflection 

■vertical force applied to tire at ground 

Tinspaning mass below Shock strut 

constants in tire force-deflection relationship 

polytropic exponent for alr-coinpreBsion process in shock 
strut 

air pressure in shock strut when fully extended 

mass density of hydraulic fluid 

air -volume of shock strut when fully extended 

shock-strut stroke 

duration of impact pulse 

angle between shock-stmt axis and -vertical 
weight of unsprung mass below shock stmt 

vertical displacement of landing-gear attachment point from 
position at initieO. contact 




m 


» 




>1 



vertical displacement of axle from position at initial contact 


Distributed Structure 
generalized coordinate for nth mode 

angle of twist of transverse station 
modal function for torsion in nth mode 

vertical dlsplacenent of elastic axis from position at 
initial contact 

modal function of elastic axis for bending in nth mode 

vertical displacement of station mass centers from position 
at initial contact 

modal function of station mass centers for coupled bending- 
torsion in nth mode 

modal amplitude of landing-gear attachment point for coupled 
bending-torsion in nth mode 

chordwise distance between elastic axis and station mass 
center 

wing span 

bending moment 

vertical component of applied landing-gear force 
natural freq^uency of first deflection mode 

polar moment of inertia of wing cross section about station 
mass center 

polar moment of Inertia of wing cross section about elastic 
axis 

radius of gyration of wing station about elastic axis 
lift force per unit length of span 
mass per unit length of sx>an 
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generalized mass for nth mode 
circular freq^uency of nth mode 
generalized force in nth mode 
shear 

natural period of nth mode 

chordwise distance between elastic axis and any arbitrary 
point 

chordwise distance between elastic axis and landing-gear 
attachment point 

spanwise distance from airplane center plane to any transverse 
station 

spanwise distance from airplane center plane to landing-gear 
station 

vertical displacement of any point from position at initial 
contact 

vertical displacement of landing-gear attachment point from 
position at initial contact 

vertical displacement of axle from position at initial contact 
virtual displacement of generalized coordinate of nth mode 
virtual wort in nth mode 


Equivalent Three-Mass System 

vertical displacement of center of gravity of spring-connected 
masses from position of initial contaet 

spring constant 

lift force acting on mass mf 

lift force acting on mass mg 
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T 


mg 

u 

Wf 


Zu 


mass acting directly on landing gear 
elastically si; 5 )ported mass 

natxaral frequency of vibration of spring-connected masses 
deflection of spring 

veight of mass acting directly on landing gear 

weight of elastically supported mass 

vertical deflection of landing-gear attachment point 

vertical deflection of elastically STJ 5 )ported mass 

vertical displacement of axle from position at initial contact 


A 

Ai 

Ol 

°Lo 

PLa 

7 

ro 

p 

Vl 

w 


Aerodynamic 

total wing area 

wing area assxmied concentrated at station i 
lift coefficient 

lift coefficient at Instant of initial contact 
lift-ciirve slope 
flight -path angle 

flight-path angle at instant of initial contact 

mass density of air 
landing speed of airplane 

total weight of airplane 
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Subscripts ; 

a aerodynamic 

f landing-gear atta chme nt, point 

g landing-gear station 

1 any spanwise station 

n pertaining to the nth. mode 

0 zero or rigid-body mode 

T at instant of initial shock-strut motion 

T at instant of maximum, landing-gear force 

max maximum 

The use of dots over symbols Indicates differentiation with respect 
to time t or t. All translations are positive downward (see figs. 1 
to 5 )* The absolute value of any term is indicated by | () |. 

ANALYSIS 


In order to study the behavior of a landing gear and a flexible 
airplane structvure as mutually interacting elements of a coupled system, 
the equations for the landing-gear force are combined with the equations 
of motion of the structure. The motions of the structure are treated by 
the mode-superposition approach, wherein the deflections of the structure 
are expanded in terms of its natural modes of vibration. The effects of 
interaction between the landing gear and the structure are Introduced by 
expressing the landing-gear force in terms of the motions of the landing- 
gear attachment point and the Artieel axle (or unsprung mass) rather than 
as an arbitrary function of time. 

Because conventional oleo-pneumatic shock struts do not begin to 
deflect until some finite time after initiad contact of the tire with 
the ground, the impact is treated in two parts, namely, the phases prior 
to and subsequent to the beginning of shock-strut deflection, where the 
initial conditions for the second phase are determined from the terminal 
conditions for the first phase. 

In the first part of the analysis, the equations for the landing- 
gear force are presented. Then, the deflections of the structure are 
expanded in terms of cotpled modes and the resulting equations of motion 
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for the system are presented in a general form. For the purpose of 
Indicating the effects of interaction, however, the system used in the 
numerical trend studies has been simplified hy restricting consideration 
of the structural deflections to the first two inodes of the expansion. 
Within the framework of this two-mode treatment, it is also shown that 
the airplane structure can he represented by an equivalent system of 
spring-connected concentrated masses, which may be used to simulate the 
effects of structural flexibility in jig drop tefets of landing gears. 


Landing-Gear Force 

An analysis of the behavior of the conventional type of oleo- 
pneumatlc landing gear was presented in reference 6 . In this study the 
mass above the landing gear was considered as a rigid body; the system 
treated therefore had two degrees of freedom and is schematically repre- 
sented in figure 1 . The analysis of the landing gear considered the 
velocity-squared damping of the metering orifice, the air-compression 
springing of the shook strut, the nonlinear force-deflection character- 
istics of the tire, and the internal shock-strut friction forces. Cal- 
culated time histories of the landing-gear forces and the motions of the 
system were in good agreement with ejqierlmental data obtained in drop 
tests. 

In the present study the rigid mass is replaced by a flexible air- 
plane structure, but the treatment of the landing gear is essentially 
the same as that in reference 6 . However, since conventional landing 
gears are inclined forward so as to minimize normal forces and bending 
moments due to the combination of vertical and drag forces, it will be 
assumed that the resultant force on the landing gear lies along the axis 
of the shock strut so that bending moments and resulting internal friction 
forces are neglected in the present analysis. 

In view of the fact that conventional oleo-pneumatic shock struts 
are preloaded with air and therefore do not begin to deflect until some 
finite time t^- after initial contact of the tire with the ground, the 

impact must be treated in two phases. In the first phase, since the 
strut is effectively rigid, the landing gear has only one degree of 
freedom and the motion of the conplete system of the landing gear and 
airplane is governed by the force between the tire and the ground. This 
ground force arises from the deflection of the tire and, in general, may 
be written as 


- 5Vg(^) (1) 

the exact variation depending on the tire force-deflection character- 
istics. Prior to the beginning of shock-strut deflection 
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Fvg = FVg(2f) (t ^ tr) (la) 

since zu = Zf . (This relationship is exact when the landing gear is 
vertical and holds very closely -vriasn the gear is inclined. ) 

The shock strut starts to deflect at the time t^- when the force 

exerted on the airplane by the shock sti*ut becomes equal to the air- 
pressure preloading force in the stmt. At this instant the free-body 
equation for the unsprung mass of the landing gear is 


“u'zfx + 


= PaoAa 


cos 9 + Wji 


(t = tr) 


( 2 ) 


Equation (2) provides the relationship between the terminal con- 
ditions for the first phase of the impact \rfiich, in conjunction with the 
solution of the equations of motion for the qomplete system prior to 
shock-strut deflection, determines the time tr when the shock-strut 

begins to deflect and, thus, the terminal values of the variables for 
the first phase of the impact, which also serve as the initial conditions 
for the second phase of the in^ct. 

After the shock strut begins to deflect, the landing gear has two 
degrees of freedom since the motions of the landing-gear attachment point 
and the motions of the unsprung mass are no longer the same. The eqm- 
tion for the vertical component of the force transmitted to the airplane 
by the landing gear after the shock strut starts to deflect is (see ref. 6) 




k __pAb |2 + p Aa 


n 


s 2 


^Cd.An^ 


^ \vo - Aas, 


cos 0 (t ^ tr) 


where 


(5) 


s = 

cos 0 

5 = , ~ 

cos 0 
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The equation of motion of the unsprung mass is 

muz’u + Fvg(zu) = IVg + ^u ("t ^ tr) (^) 


In equation (5) the first term represents the hydraulic force in 
the shock strut, where the factor s/|sj indicates the change in sign 
required between the con^ression and extension strokes. (During the 
extension stroke of the shock strut, because of the action of the rebound 
check valve or "snubber" Incorporated in most landing gears, the net 
orifice area An. will generally be smaller and the orifice discharge 
coefficient Cd will be different from the values which apply during 
the compression stroke.) The second term of equation (5) expresses the 
air-compression force in the strut, based on a polytropic pressure-volime 
relationship. In equation ( 4 ) , the force arising from the deflection of 
the tire may be expressed as Fv_(zu) = mzu^ for the usual type of pneu- 

O 

matlc tire, where m and r are constants for each regime of the tire- 
deflection process (see ref. 6 ). 


Equatio 33 s of Motion of the Airplane 

Differential equations of airplane structure .- In the mode- 
superposition approach, the structure is considered to deflect in its 
natural modes of vibration and the total displacement of any point in 
the system is the sum of the displacements of the point in all the modes 
considered. With this approach the motions are separated into functions 
which depend only on the space coordinates and functions which depend on 
the time variable. 

In the case of a impact the process is discontinuous at the 

instant tr when the shock strut begins to deflect. !Di the first phase 
of the Impact the shock strut is effectively rigid so that the motion 
of the unsprung mass of the landing gear is eBBentla,lly the same as the 
motion of the landing-gear attachment point and the force transmitted by 
the ^nnd^T^g gear to the airplane is the algebraic sijm of the ground force 
due to tire deflection, the inertia reaction of the vinsprung mass, and 
the weight of the \mspning mass. In the second phase of the impact, the 
motion of the unsprung mass is not the same as the motion of the landing- 
gear attachment point and the force applied to the airplane is governed 
by the relative motion between the landing-gear attachment point and the 
unsprung mass, as given by eq\xation ( 5 )* 

The notation employed in the analysis is Indicated in figure 2 . A 
typical transverse station located at a spanwise distance y from the 
airplane center plane is considered. The mass per unit length of span 
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is designated "by m. The translation of the elastic axis at the station 
is denoted by v} ^ is the translation of the station mass center; c is 
the chordvrlse distance between the station mass center and the elastic 
axis; and cp is the angle of twist of the station. The translation of 
an arbitrary point located at a chordwise distance x from the elastic 
axis is designated by z. The spanwise distance from the center plane 
of the airplane to the landing-gear st^ttion is indicated by yg. The 

translation of the landing-gear attachment point, or force-application 
point, is designated Zf; the distance between the landing-gear attachment 
point and the elastic axis is denoted by xj . 

In the most general case, the expansion of the deflection of the 
structiire in terms of its natural coupled modes of vibration may be 
written as 


w(y,t) = an(t)wn(y) (5) 

n=0 

and 

00 

«p(y>t) = ^ an(t)cpjy) (6) 

n=0 

whei*e the subscript n denotes the order of any mode, an is the gen- 
eralized coordinate in the nth mode, and Wn and cpn are the corre- 
sponding modal functions for bending and torsion, respectively.^ 

For later xise it is convenient to introduce expressions for the 
displacements at other points in the structure. . Since the translation 
of the station mass centers is given by ^ = w + ecp. 


00 

= 51 an(t)^n(y) 


n=0 


(7) 


where the modal function 5n = "'^n +'^n* translation of any arbi- 

trary point along the chord is given by z = w + 3tfp; therefore, 

00 

z(x,y,t) = 51 ®n('t)zn(y) (6) 

tt =0 


^he zero mode represents the translation of the airplane as a rigid 
body; therefore, wq = !• In the present analysis, rigid-body pitching 
is neglected; therefore cpo = 0. 
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■vrfaere the modal function = Wn + sapn* The translation of the landing- 
gear attachment point is given by Zf = v + XfCp; therefore. 


Zf(xf,yg,t) = ^ a^(t)ln(yg) (9) 

tt =0 

where the modal function In = Wn + xjcpn* 

By application of Lagrange’s equation and the orthogonality rela- 
tionships between coupled inodes, it can be shown (see, for example, 
refs. 8 to 10) that the equation of motion for the airplane in the nth 
mode may be written as 


Mn®n ^n^%^ “ ^n ~ 0^1^2, . . .) (lO) 


where Mu is tenned the generalized mass for the nth mode and 1® 
the generalized force, as determined from virtual-work considerations. 
For a continuous system. 



In practice the spanwlse mass distribution is often approximated 
by breaking up the distribution into discrete masses which are concen- 
trated at a finite number of stations along the span. With this approach 
equation (U) may be written as 





1 

where the subscript 


1 denotes any spanwlse station. 


(11a) 
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For the rigid-body mode (n = 0), since wq = Co “ ^ 



The relationship between and the external forces can be deter- 

mined by application of virtual-work principles. !^y definition, the 
work done in the nth mode by the generalized force acting through a 
virtual displacement of the genera.lj.zed coordinate of the mode is eq.ua! 
to the work done by the external forces acting through virtual displace- 
ment of their points of application in the mode. Thus, the virtual work 
done by the generalized force in the nth mode is 


6Wa = Qn 5Sn (l2) 

In the case of an airplane during landing the external forces are 
the distributed lift forces L(y) , the distributed weights gm(y) , and 
the force F transmitted by the landing gear. The virt^aal work done 
by these external forces in the nth mode is therefore given by 



Equating equations (12) and ( 15 ) gives the following relationship 
for Qq: 



(n = 0,1,2, . . .) (l4) 
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For the rigid-body mode (n = 0) , since ccq = 0 and Zq = ^0 = &0 = 
equation (l4) becomes 



subject to the initial conditions 


£>0(0) = 0 


and 


^(0) = 


If the airplane is assumed to be free of oscillations at the time 
of initial contact, 


®Ta(0) “ S'nC®) ~ ^ ^ 0) 


Since 


F(0) = - Wu 


equation (ill-) applied to the instant t = 0 gives 


an(0) 




(n 4 0) 


This relationship indicates that, in general, a finite static deflection 
in the flexible modes will be present at the time of initial contact. 

At any subsequent time the deflection ylll be equal to this initial static 
deflection plus an additional deflection which varies with time; 

that is, aji = an(o) + a^. • This substitution permits equation (l4) to 
be written as 
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Mnaa^. + Mno^i^an^ = - (F + Wu)ln 4 0) (^5) 

subject to the initial conditions 


In the remainder of the paper, for the sake of simplicity of notation, 
the subscript t will be dropped, with the understanding that a^^ repre- 
sents the time-varying part of the displacement of the nth mode, so that 
equation (I5) is written as 


M^an + MnC^^an — - (P + 4 


(15a) 


If the external forces are specified solely as functions of time, 
the eqiiatlons of motion for each mode of the system are uncoupled and 
can be. solved individually. However, when the external forces depend on 
the motions of the system, as in the case of the landing-gear force 
during a landing impact, the relationships between the external forces 
and the motions in the modes serve to couple the equations of motion so 
that they must be solved simultaneously. Furthermore, in the case of 
landing impact, since the process has two phases, as previously discussed, 
the equations of motion for each phase must be solved separately, where 
the initial conditions for the second phase are the same as the terminal 
conditions for the first phase. 

Motion prior to beginning of shock-strut deflection .- Since the 
shock strut is effectively rigid in the first phase of the Impact, the 
force transmitted by the landing gear to the airplane, P in eq\ia- 
tion (l5s-)j is equal to the ground force Fy (z^) less the inertia 

o 

reaction of the unsprung mass and the weight of- the unsprung mass, as 
may be seen by considering the xmspr-ung mass as a free body; thus. 


^t^tT = “u^f " 


so that the motions of the system during the first -phase of the in^iact 
are governed by the following set of differential equations: 
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Mq^ 


= - I’Vg(zf) + nHiZf - Wuj - Zq 


Mi'ai + M2CJ^^aQ_ = - Fy (z^) + 10^2^ |-j_ 

_ Q . 


1. 


^ (t ^ ti-) (16) 




FVg(2^f) + 


»m. 


where 


Zf = 

ZO = 



L dy - Mq£ 


and the mth mode is the highest mode considered. 

The initial conditions for equations (16) are the conditions at the 
instant of initial contact, namely. 


ao(0) = 0 

^(0) = Vy^ 

and 

s-nCo) = an(0) =0 (n ^ 0) 


As previously indicated, the first phase of the impact terminates 
at the time ti- when the force in the shock strut becanes equal to the 
air-pressure preload force. The terminal conditions at this instant, as 
determined by consideration of the unsprung mass as a free body, are 
given by equation (2), namely, 

m^Zf^ + Fyg^Zf^J = Pg^A^ COB 0 + 
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The solution of equations (16) in conjunction with eqmtlon ( 2 ) 
permits the detersilnation of the time tr when the shock strut begins 
to deflect and the values of the motion variables at this instant; these 
values then serve as the Initial conditions for the second phase of the 
impact. 


Motion subsequent to beginning of shock-strut deflection .- In the 
second phase of the impact the force transmitted by the landing gear, F 
in equation (15a), is the vertical conponent of the shock-strut force 

as given by equation (j) . Thus, the motions of the system during the 
second phase of the impact are governed by the following set of dif- 
ferential equations: 

MoSo = ’‘(^Yb ^ 

Miai + Miorq^aj^ = -(^^b 


= -(Fvs 


mu^ + 5 Vg(zu) = FVg + 



where 


Fvg(zf-Zu^ 


as given by equation (5); aJ^<i 


m 

^f = ®n^n 
rt=0 

The first m equations of equations (17) represent the motions of 
the airplane structure in its first m modes, whereas the last equation 
of the set is the equation of motion of the unsprung mass of the landing 
gear as previously given by equation (k -) . The initial conditions for 
equations (17) are the terminal conditions for equations (16) as previ- 
ously discussed. In view of the fact that the landing-gear forcing 
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term Py^ is highly nonlinear, analytical solution of the system, of equa- 
tions ( 17 ) does not appear possible so that it is necessary to resort to 
numerical-integration or analog methods. 


Simplified System Considered in Nimerical Studies 


The preceding section has presented the equations of motion for a 
flexible airplane coupled to a landing gear, which permit calculation of 
the motions of the system during a landing impact with consideration of 
as many modes as may be desired. For the study of the effects of inter- 
action between the landing gear and the structure, however, it appears 
that the primary effects of structural flexibility on the behavior of 
the landing gear can be represented by considering only the first deflection 
mode in addition to the rigid-body mode. 2 This simplification, which 
greatly reduces the amount of computational work, is felt to be justified 
for the puip)oses of the present Investigation since both theoretical con- 
siderations and experimental data indicate that the higher modes should 
have relatively little effect on the landing-gear performance. With this 
assumption the eq^aations of motion reduce to 

(l8a) 


( 18 b) 


Moao = - 


FVg(2f) + “u^f - Wu 


- ZO 


(t ^ tr) 




FVg(zf) + “u*Zf 


and 








(19a) 


%ai + MitD^^ai = -^Fyg + I 


(t > tr) 


( 19 b) 


™u^ + Fyg(zu) = Pyg + Wu 


(19c) 


^In a dynamic analysis, stresses in the structure due to excitation 
of the higher modes can be approximated by calculating the response of 
such modes. Individually, to the forcing function determined for the 
landing gear coupled with the rigid-body and first deflection modes. This 
procedure should be a considerable inprovement over the use of, the rigid- 
body forcing function as a basis for response calctilations in cases where 
the landing-gear attachment points experience appreciable deflections 
relative to the mass center of the system. 
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where 


Zf = &Q + 


and 


“•uZfx + = Pa^ cos e + Wu (t = tr) 

The solution of equations (18) and the determination of the condi- 
tions at the time tr when the shock strut begins to deflect, which 
serve as the initial conditions for equations (l9)# aare treated in appen- 
dix A. With these initial conditions, equations (19) may be solved by 
numerical integration or analog methods. 

From the time-history solutions for the motions of the system thus 
obtained, the accelerations and Inertia loads at any point in the struc- 
ture can be calculated from the equations presented in appendix B. 


Equivalent Three-Mass System 

It is of Interest to note that the equations of motion previously- 
presented not only represent the distributed system of the airplane but 
can also be used to define equivalent systems of spring -connected masses, 
where the number of masses above the landing gear is equal to the number 
of modes considered. For the particular case where -two modes are con- 
sidered the equivalent system is one containing three masses, one of 
which is the unsprung mass of the landing gear. The use of such a three- 
mass system provides a relati-vely simple means for simulating the primary 
effects of structural flexibility in actual drop tests of landing gears 
in a drop -test jig. 

In the equi-valent three -mass system (see fig. 3), m^ represents 
the mass to which the landing gear is directly attached and 1% is the 
elastically connected mass. The displacement of mf relati-ve to its 

position at the Instant of initial contact is denoted by Zf} the dis- 
placement of mg is designated Zg, whereas the displacement of the 
axle or unsprung mass m^ is Zg^. The spring constant of the elastic 
member is denoted by k. Separate lift forces Lg and Lf -will be 
considered to act on the masses mg and m^. 

In order that the three-mass system represents the airplane properly, 
™u^ of course, the landing -gear characteristics must be the 
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same for the two systems so that the landing -gear force is the same; and 
“■f^ applied lift forces nnist be determined ftom the 

relationships between the equations of motion for the three-mass system 
and the equations of motion for the airplane. 

Consideration of the forces acting on each mass as a free body (see 
fig* 5) leads to the following equations of motion for the three-mass 
system: 

Prior to beginning of shock-strut deflection: 


{mf + mu)zf - k(zg - Zf) + - (Wf + W^) = -IV (zf) 


mgZs + + mu)zf + (Lq + Lf) - (Wg + = -Fy (zf) 

O 


(t ^ tr) 


(20a) 


(20b) 


where 


“uZfT + ^g(zfr) = Pao^ cos 0 + Wn (t * tr) 


Subsequent to beginning of shock-strut. deflection: 


mfZf - k(zg - Zf) + Lf - Wf = -P7 


mfZf + msz's + (Lg .+ Lf) - (Wg + ¥f) = -Fy f (t > tx) (21b) 


mu.zu + Fy^(z^) = Fy^ + 


(21a) 


(21c) 


The problem is to determine the relationships between ms., mf, k, 
Lsj and Lf for the airplane so that equations (20) are equivalent to 

equations ( 18 ) and equations ( 21 ) are equivalent to equations ( 19 ) with 
the requirement that the motions of the landing gear in the three-mass 
system be the same as for the airplane, that is. 


Zf = Sq + s-qil 

and that Zu be the same in both systems. Since equations (l9c) and (21c) 
are identical, they need not be considered further in evalxiating the 
unknown constants for the three-mass system. 
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It is apparent that equations (20a) and (20h) as well as equations (21a) 
and (21b) can be written as 

m^Zf - k(zg - Zf) + Lf - Wf = -F (22a) 

mfZf + jUsZs + (Lg + Lf) - (Wg + Wf) = -F (22b) 

where 

F = FvgC^u) + “uZu - Wu 

and 


Zu = Zf 

(t % t-i-) 

Hi 

(t > tr) 


Similarly, eqtxations (l8a) and (l8b) and equations (l9a) and (l9b) 
can be written as 


+ Zq = -F 


( 23 a) 


^ ^ <c^2q^ 1 + = -F (25b) 

51 51 

Thus, the problem is reduced to determining the constants for the 
three-mass system so as to make equations ( 22 ) identically equivalent 
to equations (23). This may be done in any of several different ways. 

For example, since the structure is taken as linear, let 

Zg = 8q + a^^p 

where p is a constant to be determined. Substituting for Zf and Zg 
in equations ( 22 a) and ( 22 b) and eliminating Sq between these equations 
gives 
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“f(5i - P)ai + - P)(inf + nis)^l + Lf - Wf + -^Wg - Lg) = -F (2^1-) 

ms ms 

■whereas, subtracting equation (22a) from equation (22b), ■with the same 
substitutions, gives 

ao + pai + ^(p - 5i)ai + ^ - g = 0 (25) 

Equation (24) is directly comparable ■with equation (25h)« Combining 
equations (23a) and (23Jb) so as to eliminate F and to mate the coeffi- 
cient of So equal to unity gives the foUo^wing equation with which 
equation ( 25 ) may be directly compared: 


Ml .. Mioii^ ZO - % ^ 

ao - an - ' ■ ■ ■ ■ an + ' ■ = 0 

Moll Moll Mo 


( 26 ) 


In order to evaluate the constants for the three-mass system, each 
term in equations (24) and ( 25 ) is set equal to the corresponding term 
in equations ( 23 b) and ( 26 ), respectively, the constants in each equa- 
tion being considered as a single term. This procedure gives six simul- 
taneous equations, the solution of ■vdilch yields the foUo^wing expressions 
for the constants in "bhe three-mass sys'fcem: 


mg + ni^ = Mq 

( 27 ) 

Lg + Lf = Mog + Zq 

( 28 ) 


( 29 ) 

Ml + MqIi^ 



mf 


MiMq 

Ml + l^li^ 


( 30 ) 
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% _ 


(51) 


k 


= 


^^1 


iMi + Moll" 


(52) 


MoIi^(zq - Wu) ^ 
Ml + MqIi^ 


( 55 ) 


and 


«A + . ., 

!£. = + Wj 

Ml + MqIi^ 


where 


% 

Moll 




( 55 ) 


and 


mfZf + msZg = MoBo 


(56) 


With the foregoing substitutions, eqTiations ( 22 ) are identically 
equivalent to equations (23) j thus, the three -mass system with the speci- 
fied values of mg, mf, k, Lg, and Lf can be considered to be equiv- 
alent to the airplane in its first two modes during both the first and 
second stages of the impact. Equations (27) and (28) are required to 
satisfy the equations of motion for the airplane as a rigid body, whereas 
equations (29) to ( 3 ^) are required for proper representation of the air- 
plane in its first flexible mode. With this approach the structural prop- 
erties of the airplane are defined by three parameters: the total mass 

above the landing gear Mq, the mass ratio mg^mf, and the natural fre- 
quency £JDl. 
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The solution of the equations of motion during the first phase of 
the impact and the determination of the conditions at the instant of 
initial shock-strut deflection tr axe treated in appendix A. With 
these conditions as initial conditions, the equations of motion for the 
sepond phase of the impact can he solved by numerical- integration or 
analog methods. From the time-history solutions for the motion of the 
three -mass system, the inertia loads and bending moments at any point in 
the airplane structure can be calculated by use of the equations in 
appendix B. 


Solution of EqTiations of Motion 

In view of the fact that the equations of motion subsequent to 
time tx are highly nonlinear and therefore cannot be solved in closed 
form, it is necessary to resort to numerical-integration or analog methods. 
Various numerical-integration procedinres are given in references 11 to 13- 
Appendix A of reference 6 illustrates the application of several sxxli 
methods to the problem of the Impact of a landing gear attached to a rigid 
mass. One of these methods, which may be termed the "quadratic procedure," 
was used to obtain those numerical resialts presented in this paper which 
could not be obtained analytically. 

In this procedure, which involves a step-by-step solution of the 
equations of motion, the following difference equations (ref. 11, p. I 6 ) 
based on a quadratic variation of displacement over successive finite 
time intervals are used to replace the derivatives in the equations of 
motion: 


Tvp - 


^p+1 ~ ^p-1 

2 At 


and 


Ap+i - 2Ap + Ap_i 

(Ab)^ 

where Ap is the value of any variable at the end of the pth Interval 

subsequent to the beginning of the process and At is the time interval. 
The difference eqioations of motion obtained by substituting these expres- 
sions into the differential equations of the system then become essentially 
extrapolation formulas which peimit calculation of the displacements to 
come from the values of displacement already calculated, the whole pro- 
cedure starting out with the initial conditions of the process. With 
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the displacement time histories thus calculated, the velocities and 
accelerations are then determined from the foregoing difference equations. 


CALCULATED EESULTS AND DISCUSSION 
Cases Considered 


In order to investigate the effects of structural flexibility on 
the behavior of the landing gear and the loads in the airframe, several 
case-history studies are presented which cover a range of airplane mass 
ratios mg^mf. The calculations are based on the structural properties 
of two large airplanes having considerably different mass and flexibility 
characteristics. Airplane A is representative of a four-engine propeller- 
driven World War H bomber having a gross weight of 47,200 pounds and 
a natural frequency of vibration in the first coupled bending-torsion 
mode of 3*37 cycles per second. The structural characteristics oised for 
airplane B are representative of a present-day swept-wing six- jet-engine 
bomber having a gross weight of 125,000 pounds and a natural frequency 
of 1.29 cycles per second in the first coupled bending-torsion mode. The 
landing-gear characteristics used for airplane A were based on the man- 
ufactiarer’s data, whereas, for airplane B because information was not 
available, the shock-stmxt characteristics were chosen so as to yield a 
landing gear which is essentially a scaled-up model of the landing gear 
of airplane A. The pertinent numerical data for airplanes A and B are 
given in tables I and II, respectively; the modal functions for the first 
coupled bending- torsion mode are plotted in figure 4. 

The main landing gears of airplane A were located in the inboard 
engine nacelles very close to the nodes of the first coupled bending- 
torsion mode; in the case of airplane B the landing gear is of the bicycle 
type and is located in the airplane center plane. The position of the 
landing gear (since it determines the value of the modal amplitude Sq) 

in conjunction with the valties of Mq and Mi governs the value of the 
mass ratio nis/“if each case. (See eq. (31)-) 

In order to represent a broader range of mass and flexibility effects, 
the calculations for each airplane were made for four mass ratios corre- 
sponding to three arbitrary landing-gear positions in addition to the 
original landing -gear location. In practice, of course, a change in 
landing-gear location would, probably necessitate a modification of the 
wing structure and result in some change in the modal characteristics 
and, thus, the mass ratio. The main purpose of the calculations, how- 
ever, is to indicate the effect of mass ratio on the behavior of the 
system, and the exact locations of the lending gear which correspond to 
the mass ratios used are of secondary interest. 
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In the calcTilation of the mass ratio the landing-gear force 

was assumed to pass through the mass center of the landing-gear station. 
Since the modal characteristics used were for the complete airplane 
incli;iding the unsprung mass of the landing gear it was asstnned that 

the unsprung mass was rigidly connected to the mass mf in the equivalent 
three-mass system, as in the first phase of the impact, so that 


ms ^ 

mf + mjj 

where Mq, M^, and include the effects of the unsprung mass as 

part of the airplane mass distribution. The mass ratios considered and 
the corresponding landing -gear locations are as follows: 


Airplane A 

Airplane B 

Landing-gear 

Mass ratio. 

Landing-gear 

Mass ratio. 

location at - 

mfl/mf 

location at - 

mg/mf 

Station 0 

0 . 2 k - 

Station 0 

0.22 

Nodes 

0 

Nodes 

0 

Station 2h^ 

•52 

Station k2i0 

.83 

Station 307 

3.55 

Station yoh 

2 .&kr 


When the landing gear is located at the node of the first flexible 
mode, this mode, of course, is not excited and, since higher modes are 
not considered in the numerical calctOations , the airplane behaves as 
though it were a rigid body, its motion being governed by equation ( 23 a) . 

As might be espected, the farther away the landing gear is from the nodes, 
the larger is the effective flexibility of the system and, thus, the mass 
ratio . 

In the calctilation of the time histories of the motions of the system, 
the lift force was ass-umed to be constant during the impact and equal to 
the total weight of the airplane, that is, 

pb/2 

/ L dy = MQg + Wii 
^ 0 
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so that 


Zo = 

This assumption corresponds to the condition that 

Ls = Ws 


and 


Lf = Wf + 

in the equivalent three -mass system. 

On the basis of the calculations in reference 6, the shock-strut 
orifice discharge coefficient -was assumed as O.9 and the polytropic 

exponent n for the air-compression process vas taken as 1.12. 


Effect of Interaction on Behavior of System 

Time-history solutions for the motions of the system during impact 
at an initial vertical velocity of 10 feet per second have been made for 
the eight configurations previously mentioned. Figures 3 to 8 show the 
variation during Impact of the more important quantities, such as the 
landing-gear force F, the responses ao^g, aq/g, zf/g, Zs/s^ the 
landing-gear-motion variables, and the accelerations at the mass centers 
of several stations along the span. Comparison of the calculated results 
for the flexible cases with those for the airplane as a rigid body (or 
landing gear at nodes, mg^mf = O) indicates that the interaction between 
the flexible structure and the landing gear can result in an appreciable 
reduction in the applied landing-gear force (and thus, the nodal accelera- 
tion), the largest reductions occiirring at the highest mass ratios. Fur- 
thermore, the reductions in landing-gear force at the higher mass ratios 
are greater for airplane B, because of its lower natural frequency, than 
for airplane A. 

Consideration of the calculated time histories of the motion of the 
landing gear Indicates how the Interaction between the flexilile structvire 
and the landing, gear affects the loads produced in the' landing gear. 
Because of the flexibility of the structure, the landing -gear attachment 
point deflects upward relative to the nodes, or instantaneous center of 
mass of the system, as the applied force builds up and the deceleration 
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of the landing-gear attachment point is greater than in the case of the 
rigid airplane. Thus, the downward “Telocity of the shock-strut outer 
cylinder is more rapidly dissipated and the displacement of the outer 
cylinder is smaller throughout most of the impact. The tire deflection 
is also smaller; however, hecause of the high stiffness of the tire, the 
decrease in tire deflection is smaller than the decrease in outer-cylinder 
displacement. The net result is a reduction in strut stroke diurlng that 
part of the Impact when the maximum force occurs and an accompanying 
reduction in the strut telescoping velocity. Since the maytTmim landing- 
gear force is primarily due to the hydraulic resistance in the strut, 
heca\ise the strut stroke, and thus the air-coii 5 )ression force, is generally 
small at the time of maximum telescoping velocity, the decrease in tel- 
escoping velocity results in a decrease in shock-strut force. 

In the case of airplane A with landing gear at station 307j the 
effect of interaction is a marked change in the shape, as well as in the 
magnitude, of the time histories. Because of the superimposed vibrations 
of the structure, the shock-strut telescoping velocity (see fig. 5 ) has 
acq^tiired an oscillatory character with two peaks of the same amplitude. 
However, since the second telescoping-velocity peak occurs when the stroke 
is large, the superposition of the high air-coii 5 >resslon force on the 
hydraulic-force results in a total-force time history the second peak of 
which is much higher than the first (see F-t curves, fig. 5) and which is 
also hi^er than might he expected from the results for the smaller mass 
ratios, which have a considerably different appearance. In the case of 
airplane B, hecause of the lower natural frequency, this douhle-peaked 
characteristic does not appear even for the largest mass ratio, an mass 
ratios yielding time histories slin.i3.ar in shape, the ma.x1mum force 
decreasing in a regular manner with increasing mass ratio. 

The extent to which the first flexible modes of airplanes A and B 
are excited hy the Impacts may he observed by examining the time histories 
of a^, aj_, and aj_. As may he expected, the higher the mass ratio, the 

greater is the degree of excitation. 

From the calculated valxies of a^/g and a^/g or Zfjg and 
the acceleration at any point along the span may he computed by means of 
the equations in appendix B. Figures 6 and 8 show time histories of the 
acceleration at the mass centers of several stations for each of the 
3nndlng-gear locations considered. Becaiise of the combined effects of 
the changes in the landing-gear forcing function and in the degree of 
excitation of the flexible modes, a given change in landing -gear location 
may result in an increase in acceleration at some stations and a reduc- 
tion in acceleration at other stations. 

Figures 5 and 7 also show time histories of the acceleration Zg/g 
which would he experienced hy the elastlcaiJiy connected mass mg in the 
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equivalent three-mass system, as in a drop test. Th.e reduction in 
acceleration with increasing mass ratio is evident. As previoiisly indi- 
cated, if such a drop test were made, the measured accelerations Zf^g 

and zs/g could be used to calculate the accelerations and stresses 
that would result at any point in the corresponding airplane structiore 
by means of the equations presented in appendix B. 

Figure 9(a) presents a summary graph showing the effects of struc- 
tural flexibility and Interaction on the maximum landing-gear force for 
the various configurations considered. As previously indicated, the 
reductions in landing-gear force are greater for airplane B than for 
airplane A because of the lower natiiral frequency of the former airplane. 
For the range of mass ratios representative of existing and proposed 
large airplanes, for example, values up to about O. 3 ., reductions in 
landing-gear force up to between I 5 and 20 percent may be possible. 

Along the same lines, figure 9(b) shows the effects of Interaction on 
the acceleration response of the landing-gear attachment point and on 
the acceleration of the elastically connected mass in the equivalent 
three-mass system. 


Effects of Neglecting Interaction in the Calculation 
of Dynamic Iioads 

In the usml procedures of dynamic analysis of landing loads it is 
cvistomary to neglect the effects of interaction on the landing-gear 
forcing function and to detennine the dynamic loads in the structure by 
calculating the response of the structvire to the forcing function which 
would be obtained if the airplane were a rigid body, this rigid-body 
forcing function being either calculated or, more frequently, determined 
on the basis of drop tests of the landing gear with a rigid mass. In 
practice, either the actual rigid-body forcing function or some simplified 
analytical approximation of it (see, for example, fig. 10) is tised. 

In order to evaluate the errors introduced by neglect of interaction 
effects, the root bending moments and shears determined from the inter- 
action solutions for airplanes A and B are compared in figures 11 and 12 
with those determined by calculating the response of the various configura- 
tions to the rigid-body forcing functions previously presented and to 
simple analytical approximations to the rigid-body forcing functions. 

These bending moments and shears are total veG.ues due to both inertia 
and aerodynamic forces, the latter being included to peimit comparison 
with the steady-flight values. For reference purposes, figures 11 and 12 
also show the root bending moments and shears which would be experienced 
by a completely rigid airplane. 

The calculation of the response of systems with two degrees of free- 
dom to predescribed forcing functions is treated in appendix C. The 
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response of "bhe various configurations to the rigid-body forcing function 
was calculated by application of the numerical- integration procedure pre- 
viously described, whereas the response to the analytical forcing functions 
was obtained in closed form. The rigid-body forcing functions for air- 
planes A and B and their approximations are shown in figure 10. In the 
case of airplane A, the rigid-body forcing function was approximated by 
a pulse composed of sine and cosine segments; for airplane B, a simple 
sine pulse was used. The equations for calculating the inertia moments 
and shears from the response of the system are given in appendix C; 
simplified expressions for calculating the moments and shears d\ie to the 
aerodynamic forces are given in appendix D. 

From figures 11 and 12 it can be seen that the bending moments and 
shears calculated from the response to the rigid-body forcing function 
are larger than those determined from the interaction solutions, the 
differences being greater for the higher mass ratios where the effects 
of interaction result in a greater reduction in the magnitude of the 
landing-gear forcing function. From these particular examples, it appears 
that neglect of the effects of interaction on the landing-gear forcing 
function can lead to overconservatism in design not only of the laufiing 
gear but also of the structure, particularly for very flexible configura- 
tions with high mass ratios. As might be expected, there was relatively 
little difference in the l o ads calculated from the response to the analyt- 
ical approximations and from the response to the rigid-body forcing 
function. 

It is of interest to note that in the case of airplane A the loads 
calculated from the interaction solvrtions are greater than those calculated 
for the conpletely rigid airplane, whereas, for airplane B, the converse 
is true. This result for airplane B is due to two factors; (a) the 
dynamic anpliflcation factor is less than unity because of the relatively 
large natinral period of the airplane compared with the duration of the 
impact p^llse (ti/tn ^ 0-3)^ aJ3d (b) there is considerable reduction in 
the magnitude of the landing-gear force because of the effects of inter- 
action. In the case of aiiplane B, the natural period is of about the 
same duration as the impact pulse (ti/tn *=» 1.1) so that the dynamic magni- 
fication factor is considerably greater than unity and more than overcomes 
the effect of the reduction in landing -gear force. 

From the preceding results, it can be seen that the effects of struc- 
tural flexibility are twofold; namely, (a) a ctiange in the magnitude of 
the applied landing-gear force due to interaction, the amount depending 
on the natural frequency of the structure, the mass ratio mg^m^, and the 

landing -gear characteristics, and (b) either dynamic amplification or 
attenuation of the loads in the structure compared with those for a rigid 
body, depending largely on the ratio of the duration of the impact pulse 
to the natural period of the stmcture. In the particular examples con- 
sidered, the landing-gear force was reduced by the effects of interaction; 
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it is conceivable, however, that, for some combinations of landing-gear 
and airplane characteristics, perhaps when the natural period of the struc- 
ture is smaller than the duration of the impact pulse and the mass ratio 
is large, interaction may result in an increase in the maximum landing- 
gear force over that for a rigid airplane because of the superposition 
of oscillations of the landing-gear attachment point on the motions of 
the shock strut. Such an unfavorable effect of stnictural flexibility 
on the applied force was indicated for certain cases of seaplane impact 
in reference 7* 

In view of the foregoing observations it woiild appear worthwhile 
to consider the effects of interaction in dynamic analyses of landing 
loads when the landing gear is located at points in the airplane that 
experience appreciable deflections relative to the mass center of the 
system. 


CONCLUSIONS 


The effects of interaction between a landing gear and a flexible 
airplane stinicture on the behavior of the landing gear and the loads in 
the structvire have been studied by treating the equations of motion of 
the airplane and the landing gear as a coupled system. The landing gear 
is considered to have nonlinear characteristics typical of conventional 
gears, namely, velocity-squared damping, polytropic air-compression 
springing, and exponential tire force -deflection characteristics. For 
the case where only two modes of the structure are considered, an equiv- 
alent three-mass system is derived for representing the alrpl^e and 
landing-gear combination, which may be used to simulate the effects of 
structural flexibility in Jig drop tests of landing gears. 

As examples to illustrate the effects of Interaction, numerical 
calculations, based on the structural properties of two large airplanes 
having considerably different mass and flexibility characteristics, are 
presented. In order to cover a range of parameters, the landing gear 
of each airplane was ass-umed to be located at three arbitrary spanwise 
positions in addition to its original location. For the particular cases 
considered, it was found that 

1. The effects of interaction can result in appreciable reductions 
in the magnitude of the landing-gear force, particularly when the flexi- 
bility of the airplane structure is large and the mtural frequency of 
the structiure is small. 

2. Neglect of Interaction effects, that is, the use of the landing- 
gear forcing function for a rigid airplane in a dynamic analysis of a 
flexible airplane, can lead to the calculation of excessive loads in the 
airplane structure. 
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5. In the case of one of the airplanes, the structural loads cal- 
ciilated from the Interaction solutions are greater than those for a com- 
pletely rigid airplane treatment (rigid structure subjected to rigid- 
body forcing function) because of the fact that the effects of dynamic 
raagnificafcion more than overcome the reduction in landing-gear force dvie 
to interaction. In the case of the second airplane, because of the rela- 
tively large natural period of the structure in comparison with the dura- 
tion of the impact pulse, the dynamic magnification factor is appreciably 
less than unity. This effect, coupled with the reductions in landing- 
gear force due to Interaction, results in structural loads that are less 
than those for a rigid airplane. 

It thus appears desirable to consider the effects of interaction in 
dynamic analyses of landing loads for large airplanes, particularly when 
the landing-gear attachment points experience large deflections relative 
to the mass center of the airplane. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va. , May 5, 1955 * 



54 


NACA TN 5467 


APPENDIX A 

CONDITIONS AT BEGINNING OF SHOCK-STRUC MOTION 


Since the shock strut does not begin to deflect until the pre- 
loading force imposed by the Internal air pressure is overcome by the 
inertia forces, the shock strut is essentially rigid during the interval 
between the instant of initial contact with the ground and the beginning 
of shock-strut motion at some time t = t-,.. Dinring this interval, since 

the deflection of the tire is essentially the same as the displacement 
of the landing -gear attachment point, the system used in the numerical 
calculations to represent the airplane and landing-gear combination has 
only two degrees of freedom, namely, the rigid-body or zero-mode displace- 
ment and the deflection in the first flexible mode, the higher modes 
being neglected. The purpose of this appendix is to consider the motions 
of the system prior to the beginning of shock-strut deflection in order 
to detennlne the conditions which exist at the Instant the shock strut 
first begins to movej these motions then serve as the initial conditions 
for the equations of motion of the system during the main part of the 
impact. For this pur^se it may be reasonably assumed that the tire 
force-deflection relationship is linear for the relatively small range 
of deflection prior to the beginning of shock-strut motion and that, 
therefore, FVo.(zf) = m'zf. In order to avoid a step jump in the time- 

O 

history solution at the time t.j-, the constant m' shotild be determined 
so that 


m'zf^ 



(Al) 


Distributed System 

Prior to time t^- the equations of motion for the airplane and 
landing gear are given by equations (18) with initial conditions: 

Zf(0) = 8 ^( 0 ) = aq( 0 ) = 0 
Zf(0) = a<3(0) = Vvq 
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Since , equations (l8) can be written as 


Mq^ = -m'zf - m^Zf + Wu - Zo 


% .. ^291 

— ao + — r 


+ “u Zf + 




+ m Zf 


The exact solution of equations (A2) can be shown to be 


itjL^ Sin At - 


sin Bt 1 + 


- A^ 


0^2 -b2 


2/1 _ 1 


D cos At cos Bt + CD, 

a2 b2 ^ Vb2 A^j 


E - VE^ - ij-F 


E + VE^ - 4f 


_ + 1%) 


“l(Wu - Zd) 
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E = 


m 


(mi + MqIiS) 


+ c 


G 


F = 


G 


G = Mi(Mq + mii) + 


By successive differentiation of equation (A3 ) > the higher deriva- 
tives of Zjf(t) are found to he 




Vy^j^(A2 _ C)co8 At - (b 2 - C)cos Bt 


D 


zJt) = 


a2 - b2 


/flth 2 - b 2\ / to- 2 _ ^2' 

|sin Bt - 1 — |sln At 


(A4) 


B(b 2 - C)sin Bt - a(a 2 - C)sln At + 


^o> 2^2 » B^^cos Bt - (tu2 - a2)cos At 


(A3) 


a2 - b2 


Vv^[b2(b 2 - C)cos Bt - a2(a 2 - C)cos AtJ + 


D A^cdj_ 2 _ A2jsln At - B^oo^.^ - B2^ sin 


(a6) 
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At the time tr^ the equation of motion of the \ansprung mass of the 
landing gear as a free body is given by equation (2) which, with 
FVg = m’zf, may be written as 


“^u'zfT + m'zf^ = PaoAa cos e + 


(A7) 


Substituting for Zf^ and In equation (A7) gives a relation- 

ship between t^ and m* : 


a2 - b2 


% 


*2 2 ^ * 

— — ^-^m* - muA2)sln At - - — - muB^)sln Bt 

A B 


D 


0>,2 _ ^2 < i >,2 _ g 2 

(m' - mu.A2)cos At (m* - maB2)cos Bt + 


B^ 


^ = Pa Ag^ cos 0 + Wa (A8) 


Beca-use equation (a 8) is transcendental in both t^- and m’ 

(m* being involved in the constants A and B), in order to obtain an 
explicit solution for tr or m’, some approximation to the trigonometric 
terms is necessary, the 03rder of the approximation depending on the accu- 
racy required. For the determination of t^ and m' it will generally 
be sufficient to assume first-order approximations for the trigonometric 
terms where only the first terms of ti^lr series expansions are used. 

With these approximations the solution of eqmtion (a 8) for tx is 

cos 0 + W^ - 

bp = -i— 2 ; L (A9) 

m MiMqVYq 

As indicated previously, m' cannot be chosen arbitrarily but must be 
determined in accordance with equation (a 1) , which may be written as 

r-1 

m' = mzf^ 


m 
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The first-order approximation for obtained from, equation (A3)^ 

is 


With these substitutions equation (A9) may be -written as 


tr 


1 

G^Pa^Aa cos 9 + Wu - Dm^) 

^Vo 

mMqMQ 


1/r 


and the equation for m* becomes 


. l/r 

TTl — Tn 

G^Pso^a cos 9 + W^ ~ 

ill ^ iii 

MqMo 


r-1 

r 


(AID) 


(All) 


(A12) 


The first-order approximations for the derivati-ves of at time 

tx, from equations (a 4) to (a6 ) , are 

Zfx = (AI 5 ) 

Zfr = ’^Vo[c - (a 2 + B^)] -tr + D (A14) 

and 

’^fr = ’'^Vo[c ’ (a2 + b2^ + txDj^2 . (a2 + (AI5) 

With the values of tr and m* calculated from equations (All) 
and (A12) , the values of Zf^ = Zf^ = and Zf^ = z^^ can be 

calculated from eqTiatloriS (aIO), (A 13)7 and (a 14) , respectively. These 
values pro-vi.de two-thirds of the Initial conditions for the process s-ub- 
sequent to the beginning of shock-strut deflection (eq. 19) • The 
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remaining initial conditions, for example, 

obtained by manipulation of the differential equations (A2) . From equa- 
tion. (A2b) it can be seen that 






M- 

^l^ 


+ m I - 8 q 

Il2 / Ii2 




(AI 6 ) 


By differentiation. 


ao^ = 


&1^ 




'Ml 


\... . Ml .. 

+ zr - ao^ 


+ “u Zf^ + 
=1~ / \ 


2 ' ■ \ ^^2 ' I "r 1^2 


(AIT) 


■vdiere, from eq^uatlon (A2a) , 


sOt = - 


^ °Hi%t - ^u ^ 

Pa^ cos e + Zq 
Mq 


(A18) 


Differentiating equation (A2a) gives 



m*z. 


f + 




Mo 


(AI9) 


The substitution of equations (AI 8 ) and (AI 9 ) and_ the initial con- 
ditions previously determined ^z^^ , z^ , ^ into eq^ua- 

tlons (A16) to (A18) provides the remaining initial conditions for the 
second phase of the in 5 >act. 
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Equivalent Three-Mass System 

The equations of motion for the equivalent three -mass system prior 
to the time t-]- are equations (20) with initial conditions 


Zf(0) = Zg(0) = 0 


and 


Zf(0) = Zg(0) = VVq 

Since it has been shown that equations (20) are identically equiv- 
alent to equations (18) for the distributed system when the relationships 
between the constants of the two systems are as defined by equations (27) 
to (3^) , it follows that equations (AJ) to (AI5) are equally valid for 
the three -mass system when the constants are redefined in accordance 
with equations (27) to (3^) . The redefined constants, in terms of the 
properties of the three-mass system, may be written as 


where 


mjea^2(Mo 4- m^) 
Mo(mf + mu) 

”f(^u •* 

Mo(nif + “u) 


E 


m 


mf + mu 


+ C 


F = 


p 

m^'cD]_ 
Mo(nif + mu) 


Mq = mf + mg 
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and 


“1 = 


Imfms 


The equations for and m' , equations (All) and (A12) , become 

|l/r 


tr 




H(pap-^ cos e + 


m 


(A20) 


where 


and 




m 



cos 0 + 



(A21) 


The values of tr and m‘ given by these equations permit the 
calcvilatlon of and by means of equa- 

tions (AIO) , (AI 3 ) ^ and (Alll-) . The remaining Initial conditions for the 
second phase of the iurpact, Zg and its derivatives at the time t^-, 
can be obtained by maniptilatlon of the differential equations (20) . 
Solving equation (20a) for Zg at time tx gives 


Zsx 


= i (mf + mjzf^ + (k 


m’ ) ^fj. + Lf + Wf 


- 


(A22) 


Differentiating eqaaatlon (20a) and siibstltutlng Pv'-(zf) = m'zf gives 

O 




+ (k - 



(A23) 
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An expression for Zg^ can easily be obtained from eq^uation (20b) 
as follows ; 


"** "** + % + W^)J (A24) 


Equations (A22) to (A24) , in conjunction with the values of Zf^, 
and z^^ previovisly determined, supply all the initial conditions 
for the second phase of the impact of the equivalent three-mass system. 
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APPEMDIX B 

DYWAKEC LOADS IN AIRPLANE STRUCTURE 


The equations of motion of the airplane have been previously pre- 
sented in several forms so that solutions for the motions of the struc- 
ture can be obtained in tenns of the variables bq and e.-^, ag and Zf, 

or Zjf and Zg. The purpose of this section is to present equations 

from which the accelerations, bending moments, and shears at any point 
on the airplane structure can be calculated once the time-history solutions 
for the basic variables have been obtained. 


Acceleration 

At any point .- The absolute displacement at any point on the struQ- 
ture (see fig. 2 ) is 


z = w + jsp 


Since 


w = ag + ajWQ^ 


and 


9 = ai9i 

where w^^ and cp3_ are the modal functions for bending and torsion, 
re spe ct Ively , 


z = ag + ai(wi + xcpi) 


and 


z = ag + a^Cwj^ + 


(Bl) 
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Since 


an = 




the acceleration at any point may also be written as 


W-, + Mpi 

• • **./** •» \ sL J. 

Z = + (Zf - 

^1 

Since, from equation ( 36 ) 


(B2) 


ao = 


mjZf + mgZg 

“0 


the acceleration can also be written as 



m^Zf + fflgZg 


+ mg(^ - Zg) 


^1 + ”Pl 

61 


(B5) 


Along elastic axis .- At the elastic axis, the displacement is desig- 
nated w and X = 0 so that equation (B1) becomes simply 


«• •• •• 

w = aQ + s-qyiL 


Equation (B2) becomes 


•• •• I O \ ^1 

W = 8q + VZf - 


(Bi^) 


(B5) 


Equation (B3) becomes 
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Along station mass centers .- At the mass center of any station the 
displacement is designated ^ land x = 6 so that equation (B1) hecomes 


b* = "^ + '^1^1 (bt) 

where is the modal function for the station mass centers and is 

equal to W]_ + ecp^* liquation (B2) hecomes 

1= Kq+ (zf - (b8) 

^1 


liquation (B3) hecomes 



► 

mfZf + msZg + %(zf - ’4)-^ 

51 


(B9) 


Bending Moments 

Outboard of landing gear . - The heading moment at any spanwise sta- 
tion yj outboard of the landing-gear station yg is readily determined 

by Slimming up the inertia moments produced by the accelerations of the 
mass centers of all stations 1 between station yj and the tip. Thus, 


, ■tip 

“i^i(yi - 


(BIO) 


Inboard of landing gear .- The bending moment at any spanwise sta- 
tion yj inboard of the landing-gear station yg is equal to the sm 

of the inertia moments produced by the accelerations of the mass centers 
of all stations i between station y^ and the tip plus the moment 

produced by the landing-gear force. Thus, 
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tip 

(BM)y <y = ”5Z “I'^Cyi - yj) + 5’Cyg - yj) 
i=j 


(Bll) 


Inasmuch as 


F = - (Mq^ + 2q) 


equation (BU) can also be written as 




tip 

XI “i^’iCyi - yj) “ (Mq*^ + Zo)(yg - yj) 
i=J 


(B12) 


Shears 

Outboard of landing gear .- The vertical shear at any spanwise sta- 
tion yj outboard of the landing-gear station yg is simply the sum 

of the inertia reactions due to the accelerations of the mass centers 
of all stations 1 between station yj and the tip. Thus, 


tip 

= “X “l^’l 

Inboard of landing gear .- The vertical shear at any spanwise sta- 
tion yj Inboard of the landing-gear station yg is the side of the 

inertia reactions due to the accelerations of the mass centers of all 
stations 1 between station y^ and the tip plus the landing-gear 

force . Thus , 

tip .. 

Syjgyg = X + F 




(bi4) 


X “^l^i " + 2 q) 
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APPENDIX C 

RESPONSE TO GIVEN FORCING FUNCTIONS 


In this appendix equations are presented for the acceleration response 
of the airplane structure to predetermined forcing functions applied by 
the lanfHng gear. The cases considered are the arbitrary forcing function, 
the sine pulse, and a pulse made up of sine and cosine segments. For the 
particular case -where the landing-gear forcing function can be represented 
by a single sine pulse, 

F(t) = Fjoajj sin Gt 

■where fl is the circular frequency of the applied sine pulse and is 
expressed by 



where T is the time -bo reach Fmax* 

If the forcing pulse is not symmetrical in time about its maximi.im 
value, it may be represented by a combined pulse consisting of a sine 
function up -fco the time T and a cosine function subsequent -to the 
time T. This latter function may be writ-fcen as 


= Fjoax ^l"^’ t' > O 


where 


t' = t - T 

and Gq is the circular frequency of the cosine pvilse; the initial con- 
ditions are the same as the conditions at the time t = T determined 
from the response to the sine-function segment of the pulse. 
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The solutions are presented for the distributed system of the air- 
plane (sketch a) and for the eq^ul valent concentrated-mass system (sketch b) 



(a) Distributed system. (b) Concentrated-mass system. 


Distributed System 

The acceleration response of the rigid body or zero mode is immedi- 
ately evident from the equation of motion for n = 0 , namely, 


F(t) + Zo 
■- 

Mo 

The response of the deflection modes follows. 

Arbitrary forcing function .- When the landing -gear forcing function 
is predetermined and arbitrary, the equation of motion for the nth mode 
(eq. (15a)) can be written as 

Mn^n + MnC%^an = - [p(t) + In (n ^ 0 ) (Cl) 

where F(t) is an arbitrary function of time and an is the generalized 
coordinate of the nth mode. 
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The general solution of equation (Cl) may be written as 


an(t) ^ / F(T)sln o>n(t - t)<1t + ^ ^^( cos «%t - l) + 

Jo 






an(0)cos aijjt + 


^(o) 


sin ojQt 


(C2) 


The acceleration response is obtained by double differentiating equa- 
tion (C2) as follows: 


an(t) = ^ 
Mn 


F(t) 


oin f F(T)sin oc^(t - T)dr + cos cDnt 

Jq 


an( 0 )£% 2 cQs oa^^t - ^(0)<% sin q^t 


(C5) 


Equations (C2) and (C3) are general solutions to equation (Cl) and thus 
represent the response of any mode to an arbitrary forcing function F(t) . 
In the present study of landing impact, the Initial conditions are 


a'n(O) — 0 


and 


an(0) = 0 

Sine-pulse forcing function .- For the particular case where the 
forcing function is a sine pulse, the acceleration response, as deter- 
mined from equation (C3)^ is 


an(t) = Fiaax ^ 


In 


Mn 

- 01,2 


(n sin £%t - cun sin fit) - sin fit 


+ an(0)qn^ cos cBnt - ^( 0 ) 0 % sin o^t 
where, again, a^Co) = 0 and aji(o) = 0. 


(C4) 
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Half-sine— half-cosine pulse .- In this case the response up to 
time ~T Is given by equation (C4) . Siibsequent to time T the accelera- 
tion response, determined from eqiiatlon (C3), may be written as 


s-nC'fc*) - Fmax ~ 
Mn 




Ol^ - oin^ 


(cos ant' - cos flit') - cos fiit' 


+ an(0)t%2 
Mn 


cos ojnt' - a^(0)c% sin £%!* (C5) 


where 


t* = (t - T) ^ 1 

®-n(0) = 


an(0) = ^ 


Equivalent Concentrated^ass System 

The equations of motion for the concentrated-mass system subject to 
an arbitrary forcing function are (see eqs. (22)) 


mfZf - k(zg - Zf) + Lf - Wf = - F(t) 

mfZf + msz's + (Lf + Lg) - (Ws + Wf ) = - P(t) 
Introducing the new variable 


(C6) 


U = Zs - Zf 

penults the combination of equations (c6) Into a single equation in one 
variable : 
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nifi 


m,fU + k(l + + J 


(C7) 


where 


j = Lf - Wf - ^Ls - Ws) 
mg 


The solution of equation (C7), ty analogy with equation (Cl)^ can 
he written .as 


1 rt T 

u(t) = / F(r)sin o>,(t - T)dT + (l - cos o>,t) + 

nifCJ3_ Q w2 ■*- 


lllfGiq*- 


u(0)cos o>it + sin coit 


(C8) 


where 


2 V ^ 

cdtl = k 




By suhstituting u(t) for Zg - Zj in equations (c6) and combining, 
the following eq\iations for the responses Zg and 'if can he obtained: 


•i3(t) = - i 


— ^ P F(r)sin 0 ^(t - T)dT + ^^:r(l - cos Oht) + 

Jq mfonS 


u(0)cos a>Lt + sin ojj^t 


Ws - Ls 

ms 


(C9) 
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and 




F(t) + (Lf - Wf) 


+ k 




f F(T)sln o:h(t - T)dr + 
J A 


— ^ — (1 - coa o^t) + u( 0 )cos cu^t + 

nifOC^^ 


u(0) 

031 


sin 


(CIO) 


In equations (C9) and (CIO), u(0) = u(0) = 0 for the present 
application to landing impact. 

Sine-pulse forcing function . - For the case where the forcing term 
is a sine pulse, equations (C9) Eind (CIO) become 


Zg(t) = - 


niQ 


F^(g sto sin iit) ^ 

/ o <i\ n' X ' 


- cui^) 


u( 0 )cos o>it + sin o>it 


Ws - Ls 

% 


(Cll) 


and 




F, 


max: 


k(n sin (Uit - 0 ^ sin fit) 

nLfOJi^n^ - 0 ^ 2 ^'^ 


- sin nt 


- (Lf - Wf) + 


Jk(l - cos c^it) 

mjf.o:i2 


+ k 


u( 0 )cos o:it + sin o:^t 


(C12) 


where, again, u(o) = u( 0 ) = 0 . 
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Half-sine— -half -cosine pulse .- The response up to time T is given 
ty equations (Cll) and (C12). SxdDsequent to time T, the responses 
(eqs. (C9) and (ClO)) become 


2s(f) = - 



^maac( cos 

OJ^t ' - cos 

at') 

“s 

mf| 

s 

CVI 

H 

1 


+ ^1 - COS CUit') + 

mfCi:Q_2 


u(0)cos o^t' + sin o:^t' 


0^ 


Ws “ Ls 
ms 


(C13) 


and 


Zf(t') 


-l^max 


k(cos o:^t' - cos 




- cos ft^t* 


- (Lf - Wf) + 




Jk(l - cos o>^t') 

mf0^2 


+ k 


U(0)COS CC^t + Bln CBqt' 


(C14) 


where 


t' =. t - T = 0 


u(0) = Uqj 


U(0) = Uqj 
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APPENDIX D 

AERODIWAMIC AND WEIGBI MCMEJSTS AND SHEARS 


In appendix B eq^uations were presented for the bending moments and 
shears due to the coitblnation of the inertia forces arising from the 
accelerations of tbs masses distributed along the span and the landing- 
gear force. In the calculation of the total moments and shears, however, 
consideration must be given to the aerodynamic lift and weight forces. 
This appendix presents equations for estimating these aerodynamic and 
weight moments and shears which, although only first approximations, are 
considered sufficiently accurate for the purposes of the present study. 

If it is assumed that the lift coefficient is constant along the 
span and eqiial to the average lift coefficient of the wing Ct , the lift 

Op ^ 

force at any station y^ is equal to ^ % where A^ is the area 

assuxned to be concentrated at the station. 

The moment at any station due to the lift and weight forces at 

is 


tip 

- g ^ mi(yi - yj) (Dl) 

i=J 

If unsteady-state lift effects are neglected, the instantaneous lift coef- 
ficient is related to the lift coefficient at the Instant of Initial con- 
tact by the expression 


each station 1 outboard of station y^ 


(Ma)y^ = CJl I ^ Ai(yi - yj) 


tip 

1 =.} 


Cl = + ClJ 7 - 7o) 




Vl 


I 


Inasmuch as the total lift at the instant of contact is 
b/2 

L dy = Zq + Mog, 




Zq + MQg 


o p 2 
|aVl 



NACA TN 5^7 


55 


so that 


(Ma) 




Zq + MqS ^ 
A 




tip 

II Ai(y± - yj) - 

i=J 


tip 

g ™i(yi 

i=J 



Similarly, the shear at any station y^ is 


(D2) 


(Sa) 


Zq + MQg 


yj 


+ - ^vjl Yl 


tip tip 

II Ai - g ^ (D3) 

i=J i=J 
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TABLE I.- CHAEACTERiariCS OF AIRPLME A 


( a) Structure 


Data taken from ref. 8 


Station, 

in. 

mi, 

lb-sec2 

in. 

leai, 

lb-in. -sec2 

in. 


’'ll 

0 

28.5 


0 

-0.078 

0 

133 

16.3 

85,234 

-39.26 

-.031 

-.00084 

217 

5.27 

1,288 

0 

-.047 

-.0016 

307 

9.15 

61,717 

-62.19 

.164 

-.00183 

428 

.97^ 

536 

0 

• 37 it 

-.001^ 

548 

.686 

287 

0 

.670 

-.00187 

638 

.153 

34.1 

0 

.936 

-.00188 


Mo, 

Ml, 

fl, 


lb-see^ 

in. 

lb-sec^ 

in. 

cps . . 


61.033 

1.607 

3.365 


(b) Shock strut 
Manuf actiirer ' s data 


Ah, sq. ft 0.163 

Aa, Bq ft 0.214 

A^, sq ft 0.00173 

Vo, cu f t 0.2597 

Pao, Ib/sq ft 30,528 

p, slugs/cu ft 1.626 


(c) IMsprung mass 



Wu, lb 

Tires (one per landing gear) 
Tire pressure, lb/s*q in. . . 

m, Ib/ft 

r 


700 

56-lnch smooth contour 

70 

85,309 

1.22 
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TABLE II.- CHARACTERISDICS OF AIRPLAHE B 
(a) Stmcture 
[unpublished dataj 


Station, 

in. 

mi, 

lb-sec^ 

in. 

leai, 

Ib-ln. - Bec ^ 

in. 

Wli 

9 li 

0 

109-534 

4,475,280 

200.57 

-0.0585 

-0.000176 

84 

4.695 

3,046 

-4.65 

-.0579 

-.000187 

168 

4.920 

19,490 

- 24.20 

-.0550 

-.000204 

252 

22.177 

278,942 

- 101.22 

.0037 

-.000231 

556 

2.5^ 

2,161 

2.44 

.090 

-.000272 

420 

2.557 

1,988 

2.60 

.1842 

-.000322 

504 

1.773 

1,136 

.92 

• 3253 

-.000379 

588 

3.269 

2,474 

-14.79 

.4772 

-.000435 

672 

8.628 

8,439 

-26.88 

.6369 

-.000482 

756 

1.144 

500 

.60 

.8181 

-.000314 

840 

.520 

186 

5.48 

1.000 

-.000526 


P 

lb-sec*^ 

“O' -TT- 

Mi, arsss£ 

In. 

fl, cps . . 


(b) Shock strut 

[values estimated from generalized curves of ref 


. 6 ] 


slugs 

2 (CdAn)^' 

sq_ f"fc • . • 
Vq f cu f *fc • ■ • 
Pao^ Ib/sq ft . 


(c) Unsprung Mass 


[Manufacturer 


's data 


] 


Wu, It 

Tires (two per landing gear) 
m per landing gear, Ib/ft 



161.775 

6.9096 


1.29 


17,900 

0.585 

0.7095 

50,528 


. 2,500 
56 X 16 
280,180 

. . 1.21 















(a) System with two degrees 
of freedom. 


(h) Schematic representation of 
shock stmit. 


Figiare 1.- Dynamical system (rigid, airplane) considered in reference 6. 
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MCA TN 5^67 



(a) Airplane A. 



0 100 200 300 IjOO 500 600 700 fiOO 

station. In, 


(b) Airplane B. 

Figure 4.- Modal functions for bending and. torsion. 
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Figure 5-- Time histories of 
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of aiiTplane A dioring impact 





(a) Laading gear at station 0. (b) Landing gear at station 2^5 ■ (c) Landing gear at station 507 

Figure 6.- Time histories of accelerations at various stations along the span of airplane A. 


WACA TN 5^6t 






ft/seo If* ft/«eo 



Figure 7 .- Concluded 
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Figvire 9.- Effects of interaction- 
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Hass ratio, nig/m|. Mass ratio, mg/nif 

(b) Accelerations of landing-gear attachment point 
and. elastically connected, mass. 


Figure 9 = ■* Concluded^ 
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(a) Airplane A. 


(b) Airplane B. 


Figure 10 .- Rigid-body forcing functions and simple analytical approximations - 
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